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EXISTENCE AND UNIQUENESS OF WEAK SOUUTIONS OF THE 
COMPRESSIBUE SPHERICAUUY SYMMETRIC NAVIER-STOKES 

EQUATIONS 

XIANGDI HUANG 


Abstract. One of the most influential fundamental tools in harmonic analysis is Riesz 
transform. It maps W functions to functions for any p e (l,oo) which plays an 
important role in singular operators. As an application in fluid dynamics, the norm 
equivalence between ||Vm||/^ and ||divM||z;>+||curlM||z;> is well established for p e (1, oo). 
However, since Riesz operators sent bounded functions only to BMO functions, there 
is no hope to bound ||Vm||/,~ in terms of ||divM||/,«. + ||curlM||i«.. As pointed out by 
Hoff[5MM7. Math. Anal. 37(2006), No. 6, 1742-1760], this is the main obstacle to 
obtain uniqueness of weak solutions for isentropic compressible flows. 

Fortunately, based on new observations, see Lemma |2]^ we derive an exact esti¬ 
mate for ||Vm||/,~ < (2-i-l/A)||divM||i». for any N-dimensionalradially symmetric vector 
functions u. As a direct application, we give an affirmative answer to the open prob¬ 
lem of uniqueness of some weak solutions to the compressible spherically symmetric 
flows in a bounded ball. 


1 . Introduction and main results 

We are eoneerned with the isentropie system of eompressible Navier-Stokes equa¬ 
tions whieh reads as 

Pt + div(pU) = 0, 

{pU)t + div(pt/ ® U) + VP = pAU + {p + d)V(divC/), 

where t > Q,x e Q. c. = 2,3), p = p{t,x) and U = U{t,x) are the density and 

fluid veloeity respeetively, and P = P(p) is the pressure given by a state equation 

(1.2) P{p) = ap^ 

with the adiabatie eonstant y > 1 and a positive eonstant a. The shear viseosity p and 
the bulk one A are eonstants satisfying the physieal hypothesis 

N 

(1.3) p > 0, p + —A>0. 

The domain Q is a bounded ball with a radius R, namely, 

(1.4) Q = fifi = {xeR^; W<P<co}. 

We study an initial boundary value problem for (1.1) with the initial eondition 

(1.5) (p,U)iO,x) = (po,Uo){x), xeQ, 
and the boundary eondition 

(1.6) Uit,x) = 0, t>0,xeda, 
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and we are looking for the smooth spherieally symmetrie solution (p, U) of the problem 

(1.1) , (1.5),(1.6) whieh enjoys the form 

(1.7) p{t, x) = p{t, \x\), U(t, x) = u(t, kl)^. 

Then, for the initial data to be eonsistent with the form (1.7), we assume the initial data 
(po, Uo) also takes the form 

(1.8) po=po(\x\), Uo = uo(\x\)—. 

\x\ 

In this paper, we further assume the initial density is uniformly positive, that is, 

(1.9) Po = Po(W) > P > 0, xeQ. 

for a positive eonstant p. Then it is noted that as long as the elassieal solution of 

(1.1) ,(1.5),(1.6) exists the density p is positive, that is, the vaeuum never oeeurs. It is 
also noted that sinee the assumption (1.7) implies 

(1.10) f/(p x) + U{t, -x) = 0, xeQ., 

we neeessarily have U{t, 0) = 0 (also [/o(0) = 0). 

There are many results about the existenee of loeal and global strong solutions in 
time of the isentropie system of eompressible Navier-Stokes equations when the initial 
density is uniformly positive (refer to IfTlI 13111 511 16ll23] - l25[l28ll29H and their generaliza¬ 
tion [|20l - l22ll27ll to the full system ineluding the eonservation law of energy). On the 
other hand, for the initial density allowing vaeuum, the loeal well-posedness of strong 
solutions of the isentropie system was established by Kim [fTTl . For strong solutions 
with spatial symmetries, the authors in IfT^ proved the global existenee of radially 
symmetrie strong solutions of the isentropie system in an annular domain, even allow¬ 
ing vaeuum initially. However, it still remains open whether there exist global strong 
solutions whieh are spherieally symmetrie in a ball. The main diffieulties lie on the laek 
of estimates of the density and veloeity near the eenter. In the ease vaeuum appears, it 
is worth noting that Xin [|^ established a blow-up result whieh shows that if the initial 
density has a eompaet support, then any smooth solution to the Cauehy problem of the 
full system of eompressible Navier-Stokes equations without heat eonduetion blows up 
in a finite time. The same blowup phenomenon oeeurs also for the isentropie system. 
Indeed, Zhang-Fang ( OTlI . Theorem 1.8) showed that if ip, U) 6 CH[0, T]; H^) (k > 3) 
is a spherieally symmetrie solution to the Cauehy problem with the eompaet supported 
initial density, then the upper limit of T must be finite. On the other hand, it’s unelear 
whether the strong (elassieal) solutions lose their regularity in a finite time when the 
initial density is uniformly away from vaeuum. 

On the other hand, there are amount of literatures investigating the global existenee 
of weak solutions to the eompressible Navier-Stokes equations, sueh as ’’finite energy 
solutions” proposed and developed by Lions I1T91 . Hoff Q and Feireisl IH, ete. One 
remarkable result is due to Jiang-Zhang [fT4l . where they prove a global existenee for 
three-dimensional eompressible spherieally symmetrie flows with y > 1. However, 
whether their weak solution is unique remains a long-standing open problem. Mean¬ 
while, Desjardins [[3]| built a more regular weak solution for three-dimensional torus. 
Inspired by his work and a new observation in Lemma [2^ we will give some positive 
answer to the spherieally symmetrie flow in this paper. 
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In the spherical coordinates, the original system (11.11) under the assumption (1.7) 
takes the form 


( 1 . 11 ) 


^ Pt + {pU)r + {N - 1)— = 0, 
r 


ipu)t + [pu^ + Pip)^ + (A^ - 1)^^ = K^Ur + {N - 1)-j 
where k = 2p + A. Now, we consider the following Lagrangian transformation: 

(1.12) t = t, y= f p(t, s) s^~^ds. 

Jo 


yt = -pur^ \ r, = u, Vy = (pr^ ‘) ^ 


Then, it follows from (1.10) that 
(1.13) 

and the system (11.1 II) can be further reduced to 

Pt+ p^{r^~^u)y = 0, 


(1.14) 


r’ ^Ut + Py = K^pir^ ^u)y^ 


where t > 0, y 6 [0, Mq] and Mq is defined by 


(1.15) 


Mn 


J r*R 

Poir)r^-^dr= p{t,r) r^-^dr, 
0 Jo 


according to the conservation of mass. Note that 
(1.16) r(t,0) = 0, r{t,Mo) = R. 

We denote by Eq the initial energy 


(1.17) 


J7 - ^ ^PI\ N-lr 


and define a cuboid Qt^y, for t > 0 and y e [0, Mq], as 
(1.18) Qt,y = [0,t]x[y,Mo]. 

From now on, we denote for p 6 [1, oo) and a radially symmetric function / = /(|y|) 

1 

/ 

ll/llz/(n) 


= ( f fdxy =coAfl f f‘’r^ ^dr 

Ja \Jo 


(1.19) 
where 

(1.20) com = A^|5il, |5i| stands for the volume of N-dimensional unit ball. 

We first prove a local existence of weak solutions to the compressible Navier-Stokes 
equations in Theorem ll.il 

The initial data are supposed to satisfy (11.8111.91) and 


( 1 . 21 ) 


>0 e L“(f2) 
Uo e H\af. 
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Theorem 1.1. Assume Q. = is a bounded ball in for N = 2,3 and y > 1, then 
there exists Tq g (0, +oo] and a weak solution {p, U) to the Navier-Stokes equations 
di.il) in (0, Tq) such that for all T <Tq, 

'peL‘-((0,T)xQ), 

pU = p{Ut + U ■ VU), VU G l 2 (( 0 , T) X 

( 1 . 22 ) , 

divU eL^{0,T;L°°ia)), 

VU G L“(0, T; L^(Q)f n L^(0, T; L°°(Q.)) 

Remark 1.1. The key idea to establish local existence of weak solution with regularity 
di.22D is to derive uniform upper bound of the density, which is analogous to Des¬ 
jardins /IH/. However, there are two obstacles in bounded domain. First of all, due 
to the lack of commutator estimates, whether there is local weak solution with higher 
regularity di.22D for system di.il) remains unknown for initial boundary value prob¬ 
lem. We rewrite it in spherically coordinate and derive some new estimates to play a 
critical role as commutator estimates frequently used by / 12 ] U9^ for Cauchy problem 
and torus. On the other hand, general global finite energy weak solution with y > 1 
was proved by Jiang-Zhang /fi?!/ . However, their weak solution is much less regular 
than di.22D and leave a challenging problem on uniqueness of such weak solutions. 
The main value in Theorem li.il is to weak the assumption from y > 3 in /15|/ to y > 1 
for three-dimensional spherically symmetric flow, and further more, give an affirma¬ 
tive answer on uniqueness of weak solutions stemming from regularity class di.22l) . 
which is our main issue in Theorem \1.2\ The technical part lies on the combination 
of Caffarelli-Kohn-Nirenberg ^ inequalities with weights and pointwise estimates for 
radially symmetric functions, see Lemmas \2.1^2.2\ 

Remark 1.2. From many early works on the blowup criterion inTMT2\\2^ of strong 
solutions to the compressible Navier-Stokes equations, the uniform bound of the den¬ 
sity induces the regularity of pU, VU, divU and curlU, as indicated by the first three 
lines in Theorem li.il Besides, one of the most important observations is Lemma 12.21 
which gives desired bound for ||Vt/||z,oc in terms of\\divU\\i,‘>-. This is a key ingredient 
in proving the uniqueness of weak solutions illustrated in the following theorem. 

Theorem 1.2. Let (p‘, U‘) be two weak solutions of( 1) obtained by Theorem \Ll\ Then 
(1.23) p^ = p2, = U^ a.e on (0, T) x Q. 

Remark 1.3. There are many weak-strong uniqueness results concerning 

compressible Navier-Stokes equations. However, weak-weak uniqueness incorporates 
more difficulties and need special attention. Up to now, the most far reaching result 
appears in m. As pointed out by Hoff the main obstacle to prevent us from es¬ 
tablishing uniqueness is whether we can prove VU G L*L°° instead ofVU & Lf BMO. 
Such a fact was first verified by Hoff f^, where the Lipschitz regularity of the velocity 
U was established with piecewise C" density. This is a only result concerning unique¬ 
ness of weak solutions. With the help of Lemma \2.2\ and Theorem \Ll\ we also give an 
affirmative answer to the spherically symmetric case. 

2. Proof of Theorem 11.11 

First we recall the following famous Caffarelli-Kohn-Nirenberg |l2l inequalities with 
weights. 
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Lemma 2.1 (Caffarelli-Kohn-Nirenberg). There exists a positive constant such that the 
following inequality holds for all u G C^fR'f 

\\-a 

\l‘i 


( 2 . 1 ) 


\\A''u\^r ^ C||YrDM|^||Y|^M| 


if and only if the following relations hold: 


( 2 . 2 ) 

(2.3) 
and 

(2.4) 

satisfying 

(2.5) 


1 , r 


-h — — a \ —h 

r n 


1 O' - 1 


p n 
0 < a - cr if a > 0, 


. 1 P 

+ (1 - a) - +- 


q n 


1 Of — 1 1 y 

a - cr < I if a > 0 and — I-= —h —. 

p n r n 


p,q>l, r>0, 0<a<l 

1 a - 1 1 /3 1 7 

—I-, —I—, —I— >0, 

p n q n r n 


y = aa + (1 - a)/3 

Furthermore, on any compact set in which d2.5l) and 0 < a - cr < 1 hold, the constant 
C is bounded. 


The following lemma is essential in proving Theorems II.H and 11.21 

Lemma 2.2. Assume Q is either a bounded ball Br with radius R or the whole space 
R^, then for any N-dimensional radially symmetric vector functions U(x) = u{\x\)-^^ 
for X eO., we have the following estimates 

(1) We have 

(2.6) ^WdivUh^^a) < l|Vf/||L»(n) < (2 + ^)\\divUh^f^a>. 

(2) For p 6 [1, oo), we have 

1^1 < {^f~~^r~^o)A’\\divU\\u>(ci) 

\ur\ < |l + (A^ - l)(^)^“^r“pm/j ||c(ivt/||z/(n) 

= iV|5l|. 

Proof Set r = |y|, obviously, 

(2.9) divU = Ur + {N-1)-, curlU = 0. 

r 

Denote by 

(2.10) Ur + {N-\)-=F. 

It follows from (12.101) that 

( 2 . 11 ) 


(2.7) 

where 

( 2 . 8 ) 


{r^~^u)r = r^-^F, 
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which gives 


( 2 . 12 ) r^-^u= f s'^-^Fds < mkA f s^-^ds) = ^ 

Jo Jo ^ 

Consequently, 

(2.13) l-l < i||F|L~. 

r N 

And 

(2.14) 


Ili^llL- 


\Ur\\L- < IFIIl- + (A^ - l)li-||L» < 

r N 


On the other hand, as = -, we have 


XiXj Xj^ 


(2.15) 


dx,{^ — ) = Ur—^ + m( —);c, 

I , „ , u 

-Ur + (dij — ) —. 

fZ J yA y 


XiXj 


It immediately implies 

(2.16) l|V[/||ioo(Q) < ||M,.||z,oo(f2) + 2||-||z,oo(n) < (2 + —)||i^||L”(n)- 

For p G (1, oo), r G (0, oo), one obtains 


„v- 


'm = r 5^ ^Fds 

Jo 


(2.17) 




^ds 


l-T 


^1- 


I f N -- 


Hence, 

(2.18) 

Therefore, (12.71) follows from (12.181) and (12.91) 
This finishes the proof of Lemma. 


M<(-r~F~oj/jjFllLr(a) 


□ 


We only prove the ease when N = 3 sinee the ease N = 2is even simpler. Through¬ 
out of this seetion, we assume that (p, U) with the form (1.7) is the solution to the 
initial boundary value problem (1.1),(1.5),(1.6) in [0, T] x O, and we denote by C 
generie positive eonstants only depending on the initial data and time T. 

We give a sketeh of proof of Theorem ll.il sinee the main idea ean be borrowed from 
Desjardins 0. Denote 

(2.19) 0(0 = 1 + ||p|L»(n) + + \\yU\\l,^^y 

Our main proeedure is to derive the following estimate for t < 1, 

(2.20) O < C + C exp |c exp(^ ;r(®)4(5)^^5)j ■ 

for some positive inereasing smooth funtion and integrable function d(5). The 
existence of Tq < 1 follows immediately from (12.201) . 
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Say denoting by ^{t) the right hand side of (12.201) . we eonelude that 


( 2 . 21 ) 


< cnmmx 

at 


for some smooth funetion n. Thus, we have 

du 


( 2 . 22 ) 


du ^ p 
Jo n(u) ~ Jo 


< I A(s)ds, 


so that there exists To sueh that for all T < Tq < 1, 

(2.23) O < Ct. 

We first have the following basie energy estimate. Sinee the proof is standard, we 
omit it. 


Lemma 2.3. It holds for any 0 < t <T, 

\U\^ apy 


(2.24) 


/ 


+ 


y- 1 


dx + K 


n \VU\^dxdT< f 

Jci 


Po- 


|C/o|2 , aPo 


+ 


y- 1 


dx. 


or equivalently, 


(2.25) 


/ u^ ap^ 

1 K"— 

“Jo r2 


r^-^dr + j< 


2 7 

aPo 1 ..iv-i 


ff 


u 


ut+'—\r’^-^drdT 


P^-^dr = Eo. 


Denoted by 

(2.26) G = KdivU - P 

as effeetive flux. The momentum equations (1.1)2 can be rewritten as 

(2.27) pU = VG. 


The main diffieulty lies in the bound of the density. We will work it in spherieal 
eoordinate as system (11.141) as follows. 

Lemma 2.4. There exists a smooth positive increasing function x{x) and an integrable 
function A{t) such that 
(2.28) 

Pit, y)<C exp [cm + C exp [cm + C , (6 y) 6 [0,1 ] x [0, Mo] . 

where 

( 2 - 29 ) m = \\P\\l,^^^ + 

Remark 2.1. From the work ofl^, refer to (103), we can prove 

(2.30) m < C exp |c ^ yr(0)d(^)j. 

Therefore, the left work is concentrated on Lemma 12.41 
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Proof. Step 1. In view of (11.141) . it holds 


.Pa 


(2.31) 


K{\Ogp)ty = K{ — )y = -Kipir U)y)y = -T Uf- Py 


-ir^-^u),-py-iN-l)^. 


Thus, integrating (12.311) over (0, t) x (0,y), we deduee that 

J' z) - (r^~^u)(t, z)) dz 


1 pit,y) 1 Po(y) , 

/clog = /clog —^ + 


(2.32) 


’ p(t, 0 ) 


Po(0) 


+ {p{s,0) - p(s,y))ds 


f 


(N - 1)— —dzds. 


If 


whieh is equivalent to 

p(Py) ^ Poiy) 

p(t, 0) po(0) 


exp /c 


(2.33) 


f((r' 

Jo 

exp {pis, 0 ) - p{s, 

ff 


u)(0,z)-(r^ ^u){t 


,z))^(z| 


• exp -K 


Po(j)n3^^vj/, 


(iV-1) 


j)) ds 

u^{s,z) 

7n 


dzds \. 


Po(0) 

We first deal with 'Ti and ¥ 3 . 

Step 2. It follows from Lemma [Z2l with N = 3, energy equality (12.251) and y > 1 


that 


(2.34) 


r r^^\u\dy= f p\u\dr <C\\VU\\L 2 ^^) f ps '^ds 
Jo Jo Jo 

< C||Vt/||i2(n) P^^s^ds ni 


3y+2 \ f’2' 

s ^y-^ds 


< Cr'-^\\P\\l,^JVUh2^n) 

<c + C\\P\\l,^^^ + 


/' 


<c + cm + c I iiPiii"ci?p 


where we used the following fact. 

First recall that 

(2.35) (P% + diviP^u) + ( 6 y - l)Pdivu = 0, 

one immediately has 


(2.36) ||P||‘,,„, < + C 


ff 


/' 


jP’^divujdxds < C + C I ||P||^Jj'. 


Hence, 


( 2 . 37 ) 


'Pi,'P7' < CexpfcAO + C WPWl^ds 
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Step 3. Similarly, recall Lemma [2^ and G = KdivU - P = kF - P 


r^u 


s^Fds 


(2.38) 


<Cr^( I F^s‘^dsY 

Jo 


Hence, 


(2.39) |m|Vt < f F^s^dsp < CrH f G^sUsY + Cr\\p\t^. 

Jo Jo 

Also, we will use the following CKN inequality easily from Lemma [LT] 


(2.40) 

or equivalently, 

(2.41) 


f 


G^s^ds < CIIVGII 


L2(B(0,r))ll^llL2(B(0,r)) 


dr 


Indeed, in view of (12.11) . take n = 1,7 = |,r = 6, Qr=jS = 1 > cr = p = q = 2,a 

Now we are ready to give estimates for ¥ 3 . 

(2.42) 

Jo Jo r Jo Jo ^ 

< C j ( j F^s^dsY f ps~^dsdT 
Jo Jo Jo 

<C f r^llplU r G'^/j^)Ur^||p||^ 

Jo Jo 

+ ibii»iiv(;iiJ,iivGiiy„„^,| + ik)ii7‘)rfT 

s c Ll + Ilpl|-'||VG||J,|,|„,„ + + 1) + \\p\t*')dT 

^ 0 


< c 


= c 




+1 


Jr 


<.cf( 


l + ||p5G||2,^^(^,.^^ + C(||p||i + l)||VG||^, 


+ 


'^jdr 


The estimates of Wp'^UWp then follows similarly as (103) in Desjardins [|3||, one 
concludes that 


(2.43) 


'Fs < 1 < < CexplC I ;r(W^)J^|- 


T' 
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Step 4. We can rewrite (12.331) as 


(2.44) 

p{t, y) = P(t) nuit, y) exp |-/c ^ J' p(s, y) ds 

where 


(2.45) 


and 


(2.46) 

m,y)=Poiy)^yv3 


On the other hand, it follows from (12.441) that 


(2.47) 

which implies 




ay 

};) = —p(t,y)^exp 

K 


- f Pis, 

K Jo 


}^) ds 


= ‘^int)m,y)y, 

K 


(2.48) exp(- r p(s,y)ds\ = il + — f (P(s)^(s,y)yds 

\k Jo I \ K Jo 

Next, we are in a position to estimate Pit). First, observe that 

/^Mo j.. 

(2.49) 


1/7 


r"" dy _ r 
Jo pit,y) Jo 


nIV 

r^-^dr=—. 

N 


In view of (12.441) and (12.481) . we have 


(2.50) 


pit,y) = 


Pit)m,y) 


(l + T S!,ins)Vis,y))yds) 


Ijy' 


Then, P{t) can be estimated as 


-L 


Pit) 

pit,y) 


dy 


(l + f i!,ins)Vis,y))yds) 


\ly 


Plit,y) 


dy 


< C 


(2.51) 


J nlVlQ ^ 


+c( 22 ;)'" 


/ 


sup y) 
V 2r,o 


< CMo 
+ CMo 


K 

( \ 
) 


sup'I/ it,y) 

V Qt,o 


L'if. 


Pisyds 


supU ^it,y) 

, Qt.o 


( \ ( \ 
-1 


sup 'l/(t, y) 

V Gr.o 


sup'ZY it,y) 

V Qt.o 


f 


P(s)yds 


l/y 


























SPHERICALLY SYMMETRIC FLOWS 


II 


Using (12.491) and taking y-th power on both sides of (12.511) . we have 


(2.52) 


‘ nty < c 

Eo 


sup y) 

\ Qt,o 


+ c 


\y 


suptY \t,y) 

V Qt.q 


V / rt 


sup nA{t, y) 

V Qt.o 


f 


PCsydsj. 


Therefore, by Gronwall’s inequality, we deduee from (12.521) that 


(2.53) 


p( 0 <cj— 

\Mo 


y-l 


/ \ 
sup 'U~^{t,y) 

k Qt,o 


■ exp |cr 

Finally, reealling (12.501) . we have 


sup'Y/ \t,y) 

V Qt.o 


\y 




sup Pl(t, y) 

V Qt.o 


(2.54) 


p{t,y) < P(t) 


sup y) 

V Qt.o 


Plugging (I2.37L(I2.43I) and (12.461) into (I2.53I) - (I2.54I) . thus finishes the proof of Lemma 

m □ 

Therefore, Theorem 1 1.1 l ean be done easily from the work of lO and Lemma 


3. Proof OF Theorem 1.2 


The remaining part is dedieated to the proof of Theorem 1 1.21 
Proof. Set 

U = 1, U' = u‘-. 
p' r 

We prove Theorem 1 1.21 for the system 


(3.1) 


Vt = (r^ ^u)y, 
i-v ^ 

U,+Py=K' 


where p = av 

The two weak solutions enjoy the following regularity 

(0<p<p‘ <p, 

(3.2) \ . f 

This is equivalent to say 

' 0 < V < v' < V, 

(3.3) \ _ u‘ 




r^-^uL - e L"L“. 


It cc 


Indeed, 
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Denote by 

(3.5) 

A = v' - V^, , © = m' - U^. 

We have 

(3.6) 

|A|<2F, A(0,y) = ©(0,y) = 0. 

and 

(3.7) 

A, = (r^-'©)„ 

(3.8) 

V-l 1 .V-l 2 

r'-"'©, + (p' - p\ = -^), 


Multiplying (13.71) by A and integrating over (0,y), we have 


(3.9) 




\A\^dy<e I l(r^-^0),rdy + C(e) I lAj^dy. 


/' 


Similarly, multiplying (13.81) by '0 and integrating over (0,y), we also obtain 


(3.10) 


<6 J l(r^-^0),fdy + C(e) J 


dy + C(e) I lAI'^Jy + C(y, v)||r^ ^Myll^co | |A|^Jy. 


/' 


where we have used Young’s inequality and 

(3.11) Ip' - p'l = |(v')-^ - (v2)-^| < C(y, v)|A|. 

Recalling the fact (13.31) and (11.131) that 

(3.12) + (N- l)v— G L^L”. 

Collecting (13.9113.lOl) and choosing e small enough, Gronwall’s inequality immediately 
imply 

(3.13) A = © = 0. 
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